Dephasing in Quantum Dots: Quadratic Coupling to Acoustic Phonons 
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A microscopic theory of optical transitions in quantum dots with carrier-phonon interaction is 
developed. Virtual transitions into higher confined states with acoustic phonon assistance add a 
quadratic phonon coupling to the standard linear one, thus extending the independent Boson model. 
Summing infinitely many diagrams in the cumulant, a numerically exact solution for the interband 
polarization is found. Its full time dependence and the absorption lineshape of the quantum dot are 
calculated. It is the quadratic interaction which gives rise to a temperature-dependent broadening 
of the zero-phonon line, being here calculated for the first time in a consistent scheme. 

PACS numbers: 78.67.Hc, 71.38.-k 



The electron-lattice interaction in quantum dots (QDs) 
plays a decisive role for understanding the time evo- 
lution and dephasing mechanisms after optical excita- 
tion. High-resolution photoluminescence spectra on sin- 
gle CdZnTe QDs show a relatively narrow zero-phonon 
line (ZPL) with Lorentz broadening T depending on tem- 
perature on top of a broad bandi The latter is a super- 
position of acoustic phonon satellites. Complementary 
and more detailed information on the polarization dy- 
namics and dephasing has been obtained from four-wave 
mixing measurements in an ensemble of InGaAs QDs. 2 
The polarization decays initially very fast during a few 
picoseconds (formation of the broad band) , followed by 
a much longer exponential decay (tens to hundreds of pi- 
coseconds) reflecting the ZPL broadening. This scenario 
is referred to in the literature as pure dephasing. 3 

In theory, the working horse has been the independent 
Boson model^ which allows an exact analytical solution 
for the linear optical polarization 1 ' 5 '^ and for the photon 
echo^ of QDs. This model starts with a carrier-phonon 
Hamiltonian being linear in the phonon displacement op- 
erators, and diagonal in the confined electronic states. 
The exact solution employs the cumulant method and 
satisfactorily describes the broad phonon band, i.e. the 
initial decay of the polarization. However, the ZPL shows 
no broadening at all (no longtime decay), in clear con- 
tradiction to the experimental findings. In other words, 
the phonon cloud induced by the impulsive excitation of 
the QD relaxes quickly to a finite lattice distortion (po- 
laron dressing of the electron-hole pair), but stays con- 
stant during the subsequent time evolution. To simulate 
the effect of pure dephasing within the independent Bo- 
son model, a polarization decay can be introduced by 
hand, or a phonon damping can be added phenomeno- 
logically via the Griineisen effect i 

To improve on this unsatisfying situation, we develop 
a purely microscopic approach to the dephasing in QDs 
which contains a finite broadening of the ZPL, thus go- 
ing essentially beyond the independent Boson model. We 
take into account the non-diagonal terms of the carrier- 
phonon interaction, i.e. phonon-assisted transitions of 
electrons and holes to higher QD levels. 3 As the acoustic 



phonons which couple to the QD have energies not more 
than a few meV, these transitions are of virtual charac- 
ter only, and do not change the carrier occupation. Using 
this argument we apply a unitary transformation which 
eliminates the off-diagonal part on the price of adding to 
the diagonal one a term being quadratic in the phonon 
displacement. Being perturbative, this transformation 
preserves the virtue of the independent Boson model, 
since the new Hamiltonian allows again a separation of 
the electron and phonon coordinates. 

With optical phonons, such a quadratic coupling in 
QDs via virtual sublevel transitions has been derived re- 
cently by Uskov et al? However, the optical response 
function was calculated only up to the second order terms 
in the cumulant, giving an (approximate) expression for 
the ZPL decay. If any damping of the LO phonons is ab- 
sent, a Gauss decay of the polarization was found. In con- 
trast, earlier theoretical worki2iii*i£ on impurity-phonon 
interaction resulted in an exponential decay, exp(— Ft), in 
agreement with the experiment. The quadratic coupling 
with acoustic phonons for the ZPL decay was explored 
recently by Goupalov et alM> and Hizhnyakov et ali* but 
assumed to be diagonal in phonon momentum. This is 
not correct since the QD breaks the translational sym- 
metry. All these approaches ^iiLiijiSii&ii however, make 
use of a long-time expansion and are not able to describe 
the initial rapid decay (the broad band). 

In this contribution we find an exact solution of the 
carrier-phonon problem in QDs with quadratic coupling 
to acoustic phonons. We calculate for the first time the 
full time-dependent optical polarization and the absorp- 
tion spectrum of a QD, including both the broad band 
and the finite broadening of the ZPL. 

After eliminating the inter-sublevel coupling the trans- 
formed Hamiltonian has the following form 

H = ^^ q a+a q + {hu eh + V)\l){l\, (1) 
q 

V = 5>/ q (a q + at q ) 
q 

-^Q qq '(a q + at q )(a q '+a^), (2) 
qq' 
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The electronic Hilbert space consists of the QD ground 
state |0) and the excited state |1) having just one electron 
in the lowest conduction-band level (El) and one hole in 
the uppermost valence-band level (El). The bare transi- 
tion energy is denoted by hw e h = E g + E\ + E\. The ma- 
trix clement for the deformation potential coupling with 
longitudinal acoustic phonons in the QD reads (a — e.h) 



2 PM u 2 s V 



-D n 



(ir^aWe^'.W, (4) 



with Vw(r) and D a being the confinement wave func- 
tion and the deformation potential constant, respectively. 
Pm is the mass density, u s the sound velocity, and V the 
phonon normalization volume. In writing Eq.Q we have 
neglected (i) the excitonic Coulomb correlation which is 
of minor importance for small QDs? but could be in- 
cluded easily; (ii) any difference between phonon param- 
eters in the QD and in the barrier material leading to 
confined modes, thus dealing with bulk phonons only; 
(iii) deviations of the acoustic phonon dispersion from a 
linear isotropic one, w q = u s q. 

The linear response function, i.e. the polarization after 
a delta-pulse excitation at t = 0, is given by the dipole- 
dipole correlation function P(t) = i (d^(t) d(0)) for t > 0, 



P(t) = ie-'^" ht { Texp 



i e - l ^ ht exp 



E 



-i\ n 1 '* 



_n=l 

x(TV(t 1 )...V(t n ) 



dt\ . . . dt r , 



with T being the time ordering, and using the interac- 
tion representation of operators. The expectation value 
for the electronic system is already performed, leaving 
only the average over the phonon operators to be done. 
In the second line of Eq. J5J we have used the cumulant 
expansion^ which has only connected diagrams in the ex- 
ponent. Applying Wick's theorem, the thermal phonon 
average forces the displacement operators to appear only 
pairwise in expectation values, which defines the (mo- 
mentum diagonal) phonon propagator 



£> q (i) - (-i/n){T[a^t)+a\(t)] [a q (0) + a q (0)] ) 



= (-i/h) (2V q + l)e- <w "l*l 



(6) 



as building block. The phonon occupation enters via the 
Bose function, iV q = l/[exp(?kj q /fcBT) — 1]. In the cu- 
mulant technique, connected diagrams are those which 
do not factorize. Thus, the propagators have to appear 



in different time integrations, forming a chain-like struc- 
ture. Concentrating first on the linear interaction, there 
is only one diagram which is connected [n = 2 in Eq. (J2J)]. 
In all higher orders, the phonon propagators have their 
time integrations exclusively, and the expressions do fac- 
torize. Things are quite different for the quadratic inter- 
action since now two displacement operators share the 
same time argument (but have different momentum) . To 
form a connected diagram, they have to be distributed on 
two different propagators, ending up with a closed loop. 
Finally, a combination of linear and quadratic interac- 
tions adds another connected structure, which is a chain 
with two open ends. The linear coupling can appear only 
twice, just at the ends. 

The calculation method for the diagram series devel- 
oped in this work can be applied to any well-behaved 
function Q qq ' . Nevertheless, a great numerical simpli- 
fication can be achieved if a simple factorization of the 
dependence on q and q' holds. Physically, this case may 
appear if one excited level having the smallest energy 
distance dominates the sum in Eq. (0. Concentrating 
at the moment on the next hole level E^, we have with 
A h = El - El the product form Q qq , = M^M^J & h . 
Then, summing over momentum, only three different 
propagators appear, 

D L (t) = l/hJ2\M cl \ 2 Dc l (t), 
q 

D Q (t) = -2/A h J2\M£\ 2 D cl (t), (7) 
q 

D M (t) = y /2/(hA h )Y l M*M$D^t) . 
q 

As a result we obtain the cumulant expansion 

P(t) = i expHuW + K L (t) + K Q {t) + K M {t)\ (8) 
with 

i /"* 

K Ut) = -- / dhdhD^h-h), 
* Jo 

K Q (t) = 

n=l J0 



. dtr. 



(9) 



K M (t) 



XD Q (t 1 -t2)...D Q (tn-tl), 
. oo , t 

= oE / dt ° ■ ■ ■ dt n+\D M {ta - h) 

n=l ^° 

x Dq(t\ — t 2 ) ■ ■ ■ -Dq(<„-i - t n )D M (t n - t n+ i) 



The first term in Eq. 10 , Kl , is all what remains from the 
linear coupling, since only one diagram survives as dis- 
cussed above. This reproduces exactly the result within 
the independent Boson model. The quadratic interaction 
brings in a purely quadratic term, Kq, and a mixed one, 
Km, which both contain an infinite sum over diagrams. 
Note that the factorial in Eq. JSJ is reduced to 1/n and 
1, respectively, due to a proper counting of equivalent 
diagrams. 
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For realistic coupling strengths, a calculation of the 
sum over n term by term is hopeless. We follow another 
path which allows us to sum up the infinite sum numer- 
ically exact: With the quadratic propagator £>q(£) as 
kernel, we look for the eigenvalue problem of the Fred- 
holm integral equation 



dt 2 Ag(ii - t 2 ) Uj(t 2 ) = Aj Uj(h) 



(10) 



The kernel is complex and symmetric, but not Hermitian. 
Therefore, the eigenvalues Aj(j — 1, 2, • • • ) are complex, 
but the cigenfunctions still form an orthonormal system. 
Note that both eigenvalues and eigenfunctions depend 
parametrically on the actual time t which defines the in- 
tegration limits. Expanding the multiple integrals into 
the orthonormal system, an n-fold convolution of Dg(t) 
converts into a power A™. The subsequent summation 
over n results in a logarithm for the quadratic term, and 
leads to a geometrical series for the mixed one, 



K Q (t) 
K M {t) 
where 



-i^lnOL-A,), (11) 

3 

jy D w (12 ) 

2 2^ i _ A, ' [iZ) 

3 J 

D Mj = \ dt 1 dt 2 D M (ti - t 2 ) Uj(t 2 ) . 
Jo 



Technically, we have discretized the time integrations on 
a fine grid, thus converting the Fredholm problem to a 
standard matrix eigenvalue search. 

Up to now we have considered virtual transitions into 
one excited hole level only. It can be shown, however, 
that the factorization of Q qq ' can be still employed for 
several excited levels if they refer to different spatial sym- 
metries. E.g., if the quantum dot has mirror symmetry, 
several excited odd states of different orientation have 
vanishing cross terms in the quadratic cumulant _Dq(£), 
and contribute additively to Kq(t). According to the 
same argument, £>a/(£) = and consequently ATm(£) = 0. 
In the actual calculation we assume a spherical quan- 
tum dot with parabolic confinement potential. The wave 
functions of the lowest states, -0ia(r), are isotropic Gauss 
orbitals with variance l a . Taking artificially l e = lh = I 
allows us to take into account the next three (degenerate) 
excited states for both, electron and hole simultaneously, 
without spoiling the factorization scheme. The only 
changes are the replacement D 2 / Ah — ► D 2 / A e + D^/A/, 
in the definition of Dg(i), and an additional factor of 
three in K Q (t), Eq. fTTT) . 

In the calculations, we have taken InAs parameters^ 
p M = 5.67g/cm 3 , u s = 4.6-10 3 m/s, D c = -13.6eV, 
D c — D v = — 6.5eV. The confinement variance I is the 
only parameter characterizing the QD. The matrix el- 
ement Eq.(@} behaves like (M^) 2 cx gexp[-(gr/) 2 /2], 
which gives for the different propagators a temporal de- 
cay on the scale l/u s . Consequently, hu s /l sets the scale 
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FIG. 1: Polarization amplitude calculated for an InAs QD 
with / = 3.3 nm at T = 100 K with linear and quadratic cou- 
pling to acoustic phonons. For the dashed curve, a finite decay 
1/r = 12.3 ps has been added by hand to the linear result. 
The inset shows the (dominant) imaginary part of the eigen- 
values Aj at ujj for t = 10 ps (circles) and 20 ps (diamonds) , 
superimposed to the Fourier transform Dq(uj) (solid curve). 



for the phonon energies involved, and gives roughly the 
spectral width of the phonon satellites (broad band). For 
the QDs studied in Ref. 0, we have extracted a value 
of I w 3.3 nm. In a spherical QD with parabolic con- 
finement, the distance to the next sublevel is given by 
A a = h 2 /(m a l 2 ). To come closer to the experimen- 
tal situation of a more pancake shaped quantum dot, 
we have instead taken the experimental observation of 
A e + Ah — 65 meV and distributed between electron and 
hole according to the InAs mass ratio of 10:1. 

The calculated polarization amplitude in Fig. ^under- 
goes a quick initial decay on the scale l/u s ~ 1 ps, which 
reflects the formation of the polaron cloud around the 
e-h-pair. This state, however, is not completely stable: 
Due to the quadratic coupling, there are virtual scatter- 
ing events into higher states which distort the polariza- 
tion, leading to an exponential decay (pure dephasing). 
Note that the quadratic coupling does not change much 
the initial decay, as seen from the dashed curve in Fig.^ 
where a decay is added by hand to the linear result. 

For times much larger than the initial decay (oc l/u s ) a 
more detailed analysis is possible. The linear term alone 
gives Ki,(t) — > — Sl + i&Lt, which is a reduction of the 
polarization amplitude exp(— S) (S is the Huang-Rhys 
factor) and a purely imaginary contribution linear in time 
(polaron shift Q). The quadratic term can be analyzed in 
terms of the eigenvalues Aj . For t larger than the tempo- 
ral width of Z3q(£i), the eigenfunctions are found to oscil- 
late in the middle of the interval as Uj{t\) oc exp(iwj^i), 
while Aj — L>Q(ujj) with the Fourier transform of -Dq(£i). 
The exact position of the eigenvalues (or real frequencies 
u)j) is determined by the boundary conditions at t\ = 
and t\ = t. Numerically, we found very good agreement 
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FIG. 2: Absorption spectra of an InAs QD with carrier con- 
finement length of I = 3.3 nm at different temperatures. The 
ZPL transition energy is taken as zero of energy. Inset: Calcu- 
lated broadening of the ZPL compared with the experimental 
results by Borri et.ali— (circles). 



with 



Aj(*)wDq 



t + 2r 



(13) 



The temporal shift r is practically independent on t and 
takes care of the 'border correction' from the finite inte- 
gration interval. 

In the inset of Fig.^we plot the eigenvalues at two dif- 
ferent times, adjusting r = 0.27 ps for optimal agreement 
with Dq(lu) (solid curve). Notice that for t = 20 ps the 
eigenvalues are two times denser than for t = 10 ps, as 
dictated by Eq. (|13|) . Combining Eqs. (|11I13|) we arrive at 



K Q (t) 



-(t + 2r) 



/o 2vr 
-Sq —Tt + iflQt 



In 



l-D Q (w) 



(14) 



been derived earlier in different models of the quadratic 
electron-phonon couplingii2*ii*iSiiii4 New is the inclu- 
sion of the time-independent correction oc r which makes 
the Huang-Rhys factor S = Sl + Sq a complex quantity, 
thus adding to the ZPL Lorentzian a (small) dispersive 
contribution. 

From the complex polarization follows the absorption 
spectrum easily, a(ui) = Im J" 00 dtP{t) exp(iwt) (Fig. |5J). 
At low temperatures, the broad band is distinctly asym- 
metric and gets more symmetric as T grows. This is not 
much different from the results of linear coupling. How- 
ever, in accordance with the decay seen in Fig.^ the ZPL 
acquires a Lorentz broadening with width T(T), which 
is exclusively due to the quadratic coupling. The inset 
shows the calculated temperature dependence of the de- 
phasing rate T(T) (solid curve). Starting with T(0) = 0, 
the full curve has nearly parabolic shape before tend- 
ing to a more linear temperature dependence. The com- 
parison with the experimental data (circles) by Borri et 
al£ shows that a quantitative agreement has not been 
reached. We attribute this to the restriction of virtual 
transitions into the next degenerate dot states. An ap- 
proximate estimate for the full sum over levels gives in- 
deed a substantial increase. Then, however, states in the 
wetting layer will be important as well. 

In conclusion, we have developed a microscopic theory 
of dephasing in quantum dots accounting for quadratic 
coupling between carriers and acoustic phonons. Intro- 
ducing a new numerical method allowed us to perform 
the infinite sum of diagrams in the cumulant. This exact 
solution gives the full time dependence of the optical 
polarization. The absorption spectrum includes the 
broad band and the broadening of the zero-phonon line 
on the same footing. Our calculated results for the ZPL 
decay in InGaAs quantum dots show the same trend as 
the experiment. 2 



which has now a real part in the linear time depen- 
dence. This gives the polarization dephasing or ZPL 
width in the frequency domain. These analytical expres- 
sions for dephasing rate T and frequency shift J7q have 
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